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ABSTRACT 


Directional solidification of the pseudobinary compound semiconductor 
material Pb-j_ x Sn x Te by the Bridgman crystal growth process will be studied. 
Natural convection in the molten sample will be visualized with a novel elec- 
trochemical cell technique that employs the solid electrolyte material yttria- 
stabilized zirconia. Mass transfer by both diffusion and convection will be 
measured by detecting the motion of oxygen tracer in the liquid. Additional 
applications for electrochemical cells in semiconductor crystal growth are 
suggested. Unsteady convection in the melt will also be detected by the 
appearance of temperature oscillations. The purpose of this study is to 
experimentally characterize the overstable conditions for a Pb-|_ x Sn x Te melt in 
the vertical Bridgman crystal growth technique and use a linear analysis to 
predict the onset of convection for this system. 
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Project Description 


1 . INTRODUCTION 


Elemental semiconductor materials, especially silicon, are well under- 
stood and allow relatively easy processing, and will therefore continue to be 
dominant in electronic systems in the foreseeable future. They suffer from 
certain inherent limitations, however; chiefly because their elemental nature 
does not allow much freedom of choice in constructing semiconductor devices. 
An obvious step toward greater freedom is the use of compound semiconductor 
materials, where the desired properties can be attained by fixing the composi- 
tion of the material. One such compound semiconductor is lead tin telluride 
(Pt>1_ x Sn x Te) , a direct bandgap material that is used as an infrared detector 
and emitter (laser or LED). The desirable features of Pb^_ x Sn x Te are its 
narrow bandgap, which gives a long wavelength (2.5 to 28 micrometers), and its 
wavelength tunability, which is achieved either internally by varying the 
composition (the value of x) or externally by varying the temperature, 
pressure, or applied magnetic field. ^ Both features are important in high 
resolution spectroscopy (specifically, atmospheric monitoring), which is one 
application of this compound semiconductor material. 

Pb^_ x Sn x Te crystals are commonly grown by a variation of the Bridgman 
technique in which a molten sample in a cylindrical ampoule is directionally 
solidified by moving the ampoule in its axial direction through a temperature 
gradient. Typically, the direction of motion is vertical which allows the 
choice of whether the hotter end of the temperature gradient will be at the 
top or at the bottom of the sample. The temperature gradient is produced with 
a three-zone furnace: a hot (above the melting point) zone and a cold (below 
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the melting point) zone separated by an insulated zone that makes the tempera- 
ture gradient more linear. 

It is convenient to consider Pb.j_ x Sn x Te as a pseudobinary solution of the 
compounds PbTe and SnTe. The phase diagram ^ shows that at all temperatures 
and compositions the liquid has a higher concentration of SnTe than the 
solid. During solidification the liquid will gradually become richer in SnTe 
as the PbTe freezes preferentially into the solid. The concentration profile 
in the solid after directional solidification will be somewhere between two 
limiting cases. The first case is the most common; it is when composition 
redistribution in the liquid melt, due chiefly to convection, occurs faster 
than solidification, and any concentration change at the liquid-solid inter- 
face quickly affects the entire liquid. The liquid is considered to be fully 
mixed at all times, and the resulting solid composition is expressed by the 
normal freezing equation 


C = kC Q (1 - g) k_1 (1 ) 

where C is the concentration in the solid of one component of the solution, Cq 
is the concentration in the sample before freezing begins, g is the fraction 
of the sample solidified and k is the segregation coefficient, that is, the 
ratio of the concentration in the solid to the concentration in the liquid at 
the interface. i n pb^_ x Sn x Te the segregation coefficient is nearly 
constant. The normal freezing case is shown graphically in Figure 1. The. 
other, more unusual limiting case occurs when there is no convection in the 
liquid and composition redistribution is by diffusion only. This case has 
been examined by Tiller, et al.,^ and is shown graphically in Figure 2. 
Most directionally solidified samples follow closely the first case but it 
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would be preferable if the second case could be achieved. Most of the 
resulting solid would then be of a constant composition and would therefore be 
more easily used for mass production of semiconductor devices. 

To achieve the second case, where composition distribution is by diffu- 
sion only, convection in the liquid must be suppressed, and this requires an 
understanding of the forces that drive convection. The strongest driving 
force for convection is gravity acting on density gradients in the molten 
material. Density gradients in Pb^_ x Sn x Te arise in two ways: thermally and 
solutally. Thermal density gradients occur because the density of the liquid 
decreases with increasing temperature, so the temperature gradient used in 
Bridgman crystal growth causes a density gradient. If the bottom of the 
sample is hotter than the top, gravity opposes the density gradient, the 
liquid is thermally unstable and convection occurs. Solutal density gradients 
occur because the density of the liquid decreases with increasing SnTe concen- 
tration. As the sample freezes, the liquid adjacent to the liquid-solid 
interface becomes richer in SnTe and therefore less dense. If the interface 
is at the bottom of the liquid, gravity opposes the density gradient, the 
liquid is solutally unstable and convection occurs. The directions of thermal 
and solutal instability oppose one another, so convection occurs whether the 
sample is hotter at the top or at the bottom. ^ It is hoped that both con- 
vective forces would be greatly diminished when crystal growth is performed in 
the microgravity environment of space. To this end, crystal growth experi- 
ments are being performed both on the ground and in space. Ground-based 
experiments have exhibited normal (convective) freezing, while space-based 
experiments have so far been inconclusive. 

The preceding qualitative description can be replaced with a more quanti- 
tative understanding by application of theory and by careful experimental 
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measurement. The differential equations that are used to model fluid flow are 
well known and cases of liquid convection similar to Bridgman crystal growth 
have been discussed in the literature; (6) however, analytic solutions to the 
fluid flow equations have not been possible without excessively simplifying 
assumptions. Current research includes numerical approaches to mathematical 
modelling of convection in the liquid. The accuracy of an analytic or 
numeric mathematical model is judged by comparison with experimental 
results. Much of the experimental work is limited to post-growth analysis of 
the resulting solid sample, which includes chemical etching to reveal the 
shape and position of the solid-liquid interface during growth, measurement of 
the solid composition by electron microprobe analysis, and X-ray diffraction 
studies to determine the orientation and quality of the crystals. Measurement 
of convection during growth has until now been limited to the use of thermo- 
couples to follow changes and oscillations in the liquid temperatures.^-^ 

In compound semiconductors the quality of the final material is very 
dependent on composition variations in the molten sample during crystal 
growth. It is proposed that convection during Bridgman crystal growth be 
monitored with thermocouple and electrochemical concentration cells. It is 
hoped that data from these measurements could be used to follow convection 
flow patterns in the liquid. The formulation of the analytical problem is 
first presented. The experimental progress to date and proposed work is then 
summarized. 

2. FORMULATION OF ANALYTICAL PROBLEM 


The system that will be analyzed is shown in Figure 3. The following 


assumptions will be made in this study: 
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1) the gas phase is inviseid and the liquid phase is a binary Newtonian 
fluid, 

2) viscous dissipation is negligible, 

3) density is a linear function of temperature and concentration, 

4) the surface tension is sufficiently high to assume a flat interface, 

5) morphological stabilities will not be considered, 

6) the temperature and velocity fields are axisymmetric, 

7) heat flow is 1-dimensional in the axial direction, and 

8) Marangoni convection can be neglected. 

The last three assumptions are especially important. The analysis of 
Schliiter, Lortz, and Bussed predict that 2-dimensional rolls are the only 
stable configuration for buoyancy-driven convection in an infinite layer of 
fluid with temperature independent properties. Also, Hoard, Robertson, and 
Acrivos^®) observed a symmetric roll pattern for deep layers of fluid with 
highly temperature dependent viscosity. Thus, it appears likely that flow 
will develop in an axisymmetric manner. 

One-dimensional heat flow is necessary for obtaining a possible conduc- 
tive state without radial gradients at the interface. Several steps have been 
taken to ensure axial heat flow: a large thickness of Shuttle tile insulation 

is used between the thermal zones of the furnace (Fu and Wilcox^ ]^), insula- 
tion has been added to the sidewalls of the quartz to cut down on radial heat 
transfer, and efforts have been made to extend the constant gradient zone and 
decrease the radial heat flux from the hot zone to the end of the sample. 
Lastly, Marangoni effects on Earth are generally dominated by buoyancy 
effects, except for very thin layers, and can be neglected completely for the 


hot on top case 
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The governing equations for this problem in non-dimensional form are: 


Continuity 


V«V = 0 


( 2 ) 


Motion 


1 3V 

— — + vVV = -VP 
Pr 3t 


Ra - - 2 

♦ ^ T F + t 1 C F + V V 
Pr Sc 


(3) 


Energy 


Pr 3T 
Sc 3t 


+ Pr (V«VT) = V T 


(4) 


Species 


|| + Sc (A-VC) = V 2 C 


(5) 


V, C, and T are the velocity, concentration, and temperature field. F is a 
unit force in the direction of gravity and t is time. These equations have 
been obtained by using the following characteristic variables: 

V - V'L/\> 

P - P'L 2 /(pC p L 2 ) 

c - (c'-c 2 )/(c 1 -c 2 ) 

T - (T'-T 2 )/(T 1 -T 2 ) 
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t = t'D/L 2 

where primed variables are dimensional and T-, and T 2 are temperature at the 
top and bottom surfaces, respectively. 

Care should be taken in choosing the characteristic variables. For 
instance, the Pr/Sc term that preceeds the time derivative in eqn. (4) is 
approximately 0.004 (Pr=.1, Sc=25). As a result, almost all of the time 
dependence of the temperature field is supplied by the interrelated motion (3) 
and species (5) equations. Thus, one could infer that the temperature field 
and concentration will be kept "in phase," with the fluctuations in tempera- 
ture following closely those of concentration. If the time had been scaled 
with thermal diffusity, though, the Pr/Sc term in eqn. (4) would not have 
appeared and the above relationship would have gone unnoticed. Instead, a 
Sc/Pr term would have preceeded 3C/8t with relatively little insight gained. 
The purpose of this study will be to use Green's functions in order to solve 
the problem in the form of an integral equation. 

3. EXPERIMENTAL DETERMINATION OF TEMPERATURE OSCILLATIONS 


The ampoule used in these experiments is a vertical circular cylinder of 
fused silica (quartz). Preliminary furnace experiments have been run with 
molten tin as the fluid phase in a destabilizing temperature gradient. A 
small amount of purified argon is trickled through the cell to prevent contam- 
ination of the tin with oxygen. Insulation is placed on the outside wall of 
the cell in order to reduce radial temperature gradients. Cotronics ceramic 
board is used to make insulating spacers that hold the Cotronics insulation 
blanket in place. The insulating board and blanket have thermal conductivi- 
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ties of 3.28 x 10 -i< cal/°C-cm-sec and 4.83 x 10* 1 * cal/°C-cm-sec, 
respectively. The quartz, on the other hand, has a much higher thermal con- 
ductivity of 7.75 x 10-3 cal/°C-cm-sec. Furthermore, the cell is placed in 
the Shuttle tile (k = 3.90 x 10 -1 * cal/°C-cm-sec) insulation zone of the fur- 
nace. Heat flow is made more 1 -dimensional in this zone due to the long 
insulation zone (65 mm or 8 r). Finally, isothermal heat pipes are used in the 
heating zones of the furnace to reduce temperature roll off at the Shuttle 
tile-heater block boundary. Reducing this roll off lengthens the constant 
gradient region which reduces the inward heat flux from the hot zone into the 
side of the cell (Rouzaud, Camel, and Favier).^^ The combination of these 
three effects should ensure axial heat flow through the cell. 

Chromel-alumel (type K) thermocouples are placed between the quartz wall 
and the insulation blanket at one cm intervals in the vertical direction. The 
thermocouples are used to measure temperature, the temperature gradient in the 
melt, and to detect temperature fluctuations (hence, unsteady convection). 
The temperatures are recorded at 5-minute intervals and one of the channels is 
monitored by a Data 6000 recording waveform analyzer at 2-second intervals. 
An FFT can be run on the data 

The thermal gradient is slowly increased until the temperature fluctua- 
tions appear on the Data 6000. When oscillations are detected, the Rayleigh 
number is calculated. The temperatures used in these calculations have not 
yet been corrected for the temperature drop that is caused by the quartz wall 
(probably on the order of 10°). Since the relevant physical properties of tin 
have a weak temperature dependence for T > 260 °C and the gradient measured in 
this way should be close to the value inside the tube, the value calculated 
for the Rayleigh number should be a reasonable approximation. Critical Ray- 
leigh numbers (Ra°^) calculated for AR = 3.15 and AR = 5.0 are approximately 
135,000 and 220,000 respectively. 
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Further experiments will involve taking temperature measurements from a 
Pb.gSn^Te crystal growth system to determine Ra° 2 and Ra s c2 . This data will 
be used to qualitatively verify results from the theoretical analysis. Other 
important experiments would include determining the density change of LTT upon 
freezing (e=1-p/p) and determining the extent of the radial temperature 
gradient. The importance of radial gradients has been discussed previously. 
One method of measuring 3T/3r will be to place a thermocouple in a capillary 
tube in the middle of the tin cell at the same axial position as a thermo- 
couple placed on the side of the cell. This cell will be heated from above, 
since the onset of axisymmetric convection in the thermally unstable case 
would produce a misleading radial gradient. 

If e is sufficiently large (say greater than .05), it can be an important 
cause of convection. Thus, it is necessary to measure accurately to determine 
whether or not to include this effect in future models. One proposed method 
of measuring e involves using a laser to mark the solid-liquid interface in a 
see-through furnace and recording the movement with a VCR camera. 

4. EXPERIMENTAL DETERMINATION OF CONCENTRATION OSCILLATIONS 


The electrochemical cells use a solid electrolyte technique that has been 
applied in a variety of thermodynamic studies. ^3) The most popular solid 
electrolyte has been zirconia (Zr0 2 ) stabilized by the addition of oxide 
dopants such a calcia (CaO), magnesia (MgO), or yttria (Y 2 0g). These dopants 
create anion vacancies in the zirconia lattice and at sufficiently high 
temperatures (>600°C) oxygen anions (0 ) can move through the electrolyte by 
migrating from vacancy to vacancy, while zirconia’s low porosity prevents the 
transfer of other species. Yttria-stabilized zirconia (YSZ) has been 



chosen for the proposed experiments because of its availability and its low 
electronic conductivity. ^ ^ 

The ionic conductivity of the YSZ solid electrolyte can be used in an 
electrochemical cell to measure the thermodynamic activity of oxygen in liquid 
Pt>1_ x Sn x Te. The cell is represented by 


M(s),M 0(s) | YSZ | Pb. Sn Te(1 ),[0], 

HI I A A 

reference electrode sample electrode 


( 6 ) 


where M and t^O represent a metal and its oxide for which the Gibb's free 
energy of formation, AG^°(M m O), is accurately known. The partial pressure of 
oxygen on the reference sie of the electrolyte, PCC^) (reference) , is known 
from the equilibrium relationship 


AG °(M 0) = RT In P(0_, reference) 1/2 , 
f m d 


(7) 


where R is the gas constant, T is the absolute temperature, and the partial 
pressure is measured in atmospheres. In eq. 6, [£] represents atomic oxygen 
dissolved in the liquid Pb-|_ x Sn x Te sample; its thermodynamic activity is 
a 0 (sample). The Nernst equation relates the thermodynamic properties of the 
cell to the potential, E, that can be measured across the cell: 


E 


RT 

2F 


In 


1 /2 

(P(0 2> reference) 

1/2 ’ 
a Q (sample) . 


( 8 ) 


F representing the Faraday constant. 


(16) 


Equations 7 and 8 can be combined to 


obtain 



( 9 ) 


E = ( 1 /2F ) [AG 0 (M 0) - RT In a. (sample)], 
f m 0 

Since a 0 (sample) is a function of the oxygen concentration in the liquid, the 
measured cell potential reveals the oxygen concentration. The meaning of eq. 
9 can also be reversed: an arbitrary potential applied across the electro- 

chemical cell will pump oxygen ions through the electrolyte until an equili- 
brium is reached, with the oxygen activity and concentration in the sample 
corresponding to the applied potential. If the cell and the applied potential 
are in a series circuit, the number of ions transferred through the electro- 
lyte can be found by measuring the number of electrons that travel through the 
circuit; that is, by measuring the current in the circuit . ^ 1 This is the 
electrochemical titration technique. It is proposed that oxygen be used as a 
tracer element in a Pb-|_ x Sn x Te sample; that the tracer be added to the sample 
by titration with one cell and detected by a second cell. The response at the 
detection cell to an input at the titration cell as a function of time will 
then reveal the nature of the convective flow in the sample between the two 
cells. 

The solubility of oxygen in Pb 1 _ x Sn x Te is finite and if it is exceeded 
the elements in the sample begin to form oxides. The first oxide to form is 
the one with the lowest AG f °, tin(IV) oxide (Sn0 2 ). Once the oxide is formed 
the activity and the partial pressure of oxygen in the sample electrode are 
fixed by the equilibrium of the reaction 

S "(1) * °2(g) ' Sn °2(s)' 

AG f °(Sn0 2 ) = RT In [a Sn P(0 2 , sample)]. 


( 10 ) 


Equation 4 can also be written as 


E 



1 /2 

P(0 2> reference) 
P(0 2> sample) 1/2 


( 11 ) 


and eq. 10 and 11 combined form 


E = (1/2F)[AG_°(M 0) - 1/2 AG.°(SnO„) + 1/2RT In a c ]. (12) 

I m f 2 Sn 


Thus the cell potential in this case can be related to the activity and hence 
the concentration of tin in the Pb i- x Sn x Te sample. Composition profiles are 
measured in solid samples after crystal growth; the electrochemical cell can 
provide data on the liquid composition during growth. 

An important part of the preparatory work for the proposed experiments 
has been the selection of the materials that will be used in the electro- 
chemical cells. An early choice was the chemical system to be used as the 
reference electrode. Gas electrodes (air, O2, or a gas mixture) are commonly 
used in solid electrolyte electrochemical cells, ^3,17) we re in this case 
rejected because of the difficulty involved in keeping the molten liquid 
sample electrode separate from the gas reference electrode. A solid reference 
electrode can be more easily assembled and sealed, so the initial choice for 
the reference electrode was the often-used nickel-nickel oxide (Ni-NiO) 
system. ^ 9-22) S ytem is well-defined thermodynamically and has a partial 
pressure of oxygen relatively near that expected in the sample electrode, 
which prevents excessive oxidation of the sample. ^3) However, the resistance 
of the Ni-NiO system to electric current and oxygen transfer, which leads to 
error in the cell potential measurement, is more than is desired and the 
copper-copper (I) oxide system (Cu-Cu 2 0) is being considered as an alterna- 
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tive.^ - ^ The next choice was the material for the electrical contacts 
from each electrode to the measuring instruments. It is essential that the 
contact materials and the electrodes be sufficiently inert towards each other 
for long periods of time at high temperatures. Experiments have shown that 
tungsten wire can be used in both liquid tin and in Ni-NiO, while the best 
contact material for Pb.j_ x Sn x Te appears to be graphite. Experiments with Cu- 
Cu 2 0 have not yet been conducted. The accuracy and effective lifetime of a 
cell are greatly enhanced if it is not allowed to oxidize. To this end the 
cell is sealed away from air and is ventilated with argon gas which has passed 
through a three-step purification line: first, hydrogen is catalytically 
reacted with oxygen to form water, then all water is removed by a dessicant, 
and finally the remaining oxygen is reacted with titanium at 850°C. The 
resulting argon has an oxygen partial pressure of about 10~33 atm. while the 
expected oxygen partial pressure of the electrochemical cell is about 10~^ 
atm. 

The proposed experiments require two types of electrochemical cells. 
First a single cell, shown in Figure 4, will be used to measure the electrical 
potential (and therefore the thermodynamic activity of oxygen in the liquid) 
as a function of oxygen concentration. The measurement will be done by elec- 
trochemically removing as much oxygen as possible from the liquid sample, then 
titrating oxygen back in while observing the changing cell potential. Once 
the liquid is saturated the solubility of oxygen in the liquid is known. The 
results from the single cell experiments form the data base required for the 
convection experiments using a double cell, shown in Figure 5. The idea for 
the double cell arrangement is from Otsuka and Kozuka.^®^ In the double cell 
the sample is in a vertical cylinder, with an electrolyte at the bottom face 
that functions as an oxygen pump cell and an electrolyte at the top face that 


functions as an oxygen detector cell. Transport of oxygen through the liquid 
sample between the electrolytes is measured by observing the response at the 
detector cell to a step or pulse input of oxygen at the pump electrode. The 
measurement will be made under two conditions: an isothermal condition where 
convection is minimal and oxygen transport is chiefly by diffusion, and a 
temperature gradient condition, where oxygen transport is by diffusion and 
convection. The first condition will permit the calculation of the diffusion 
coefficient of oxygen in the liquid (20) and the second condition will, by 
comparison to the first condition, reveal the bulk effect of convection in the 
liquid. 

The proposed experiments will give a thorough understanding of the acti- 
vity, solubility, and diffusivity of oxygen in Pb.j_ x Sn x Te and can form the 
basis for further research. For example, in Bridgman crystal growth the shape 
of the solid-liquid interface has a strong effect on the convective flow 
patterns in the liquid, ^ but it is difficult to discover experimentally the 
shape of the interface. By using the electrochemical cell as an oxygen pump 
during crystal growth, the liquid can be doped with oxygen just long enough 
for an oxygen-rich layer of solid to be deposited at the interface. After 
solidification is complete the oxygen in the solid can be detected and the 
shape of the oxygen-rich layer will show the shape of the solid-liquid inter- 
face. Another application of electrochemical cells is to use small (less than 
7mm diameter) zirconia tubes as point detectors of oxygen. These detectors 
are small enough to be located at various sites in the liquid to more 
thoroughly map the convective flow patterns in the sample. There is also 
potential for industrial applications of electrochemical cells, where they may 
be used for precise control of the concentration of oxygen as a dopant in 


semiconductor materials. 


It would also be very interesting to test the LTT system to see if the 
temperature and concentration are locked "in phase" as the governing equations 
suggest. For this, it will be necessary to measure concentration as well as 
temperture. Until now, this has been difficult in semiconductor melts. It is 
hoped that the new electrochemical cell technique can be used to measure the 
concentration of the melt while temperature readings are being made with 
thermocouples. 

Low gravity crystal growth also introduces new possibilities and some new 
problems. For example, containerless crystal growth is a promising method for 
space processing. Buoyancy-driven flow will be eliminated due to the low 

-■3 

values of g( — 1 0 g . ) and crystal defects due to wall-melt interactions 

cdrtfl 

will be eliminated. But the resulting gas-liquid surface, in the presence of 
a temperature gradient, will be subject to Marangoni convection. On the other 
hand, in conventional growth with side, walls, Marangoni convection will not be 
a problem, but convection may result from the density change upon freezing. 
If e is large, there will be a net flow of fluid towards the interface. This 
flow will not be uniform, due to the no-slip condition at the wall. So even 
though buoyancy is nearly eliminated there are still driving forces for con- 
vection in space. Perhaps the melt can be encapsulated in another liquid 
phase. The lower shear at the phase boundary will allow nearly even fluid 
flow to the interface and there will be no free surface available for Maran- 
goni convection. This is just one example of the type of problem that may be 
investigated in the future. 
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Figure 4. Single electrochemical cell assembly. 



Figure 5. Double electrochemical cell assembly. 
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